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We examine nonequilibrium features of collapse behavior in model polymers with competing crystallization and glass transitions 
using extensive molecular dynamics simulations. By comparing to "colloidal" systems with no covalent bonds but the same 
non-bonded interactions, we find three principal results: (i) Tangent- sphere polymers and colloids, in the equilibrium- crystallite 
phase, have nearly identical static properties when the temperature T is scaled by the crystallization temperature T cryst ; (ii) 
Qualitative features of nonequilibrium relaxation below T cryst , measured by the evolution of local structural properties (such as 
the number of contacts) toward equilibrium crystallites, are the same for polymers and colloids; and (iii) Significant quantitative 
differences in rearrangements in polymeric and colloidal crystallites, in both far-from equilibrium and near-equilibrium systems, 
can be understood in terms of chain connectivity. These results have important implications for understanding slow relaxation 
processes in collapsed polymers, partially folded, misfolded, and intrinsically disordered proteins. 
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1 Introduction 

Collapse transitions of single chain polymers induced by 
changing control parameters such as temperature or solvent 
quality yield rich nonequilibrium behavior when the rate at 
which these control parameters are changed exceeds charac- 
teristic (slow) dynamical rates. Investigating the glassy dy- 
1 namics of polymer collapse is important for understanding e.g. 
crystallization kinetics and protein misfolding, yet the ma- 
jority of studies have focused on equilibrium behavior. In 
this manuscript, we characterize the nonequilibrium and near- 
equilibrium collapse and crystallization dynamics of single 
flexible polymer chains. 

We employ a minimal model that yields competing crys- 
tallization and glass transitions. Monomers are modeled 
jas monodisperse tangent spheres with hard-core-like repul- 
sive and short-range attractive interactions. Recent studies 1,2 
have shown that in equilibrium, model polymers with narrow 
square-well interactions exhibit direct "all-or-nothing" crys- 
tallization transitions that mimic the discrete folding transition 
observed in experimental studies of proteins. 3 Short-range at- 
tractions also give rise to degenerate, competing ground states, 
which kinetically hinder collapse to equilibrium crystallites. 
The associated rugged energy landscapes are believed to con- 
trol the behavior of intrinsically disordered proteins. ~ 

A novel aspect of our work is quantitative comparison of 
polymer collapse dynamics to that of 'colloidal' systems with 
the same secondary interactions but no covalent bonds. Poly- 
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mers are distinguished from other systems by their topology ; 
connectivity and uncrossability constraints imposed by cova- 
lent backbone bonds give rise to cooperative dynamics- and 
phase transitions 8 not present in nonpolymeric materials. Our 
choice of tangent monodisperse spheres yields identical low- 
energy states for polymers and colloids, but (because of the co- 
valent backbone) very different free energy landscapes. 9 This 
greatly facilitates a robust comparison of crystallite formation 
and growth dynamics that isolates the role of topology and 
allows us to isolate and quantify the contributions of the cova- 
lent bonds and chain uncrossability to cooperative rearrange- 
ments and slow dynamics during collapse. We perform this 
comparison using extensive molecular dynamics simulations 
of thermal-quench-rate-dependent collapse and post-quench 
growth of polymeric and colloidal crystallites. 

Figure [T] depicts rate-dependent collapse behavior of sys- 
tems interacting via hard-core-like repulsions and short-range 
attractions. 11 In the limit of slow quench rates \T\ < |T*|, 
where T* is a critical quench rate, finite systems exhibit a 
first-order-like transition from a high-temperature "gas" (for 
polymers, a self-avoiding random coil) phase to crystallites. 
The equilibrium transition occurs if \T\ is small compared 
to key relaxation rates, such as the crystal nucleation rate r 
and rate s of large rearrangements in compact structures. At 
larger |T|, systems fall out of equilibrium, and pass onto the 
metastable liquid branch. If T becomes low enough such that 
s(T) ^> \T\, the systems become glassy and disorder is frozen 
in at T = T g i ass . Otherwise, systems relax towards equilib- 
rium crystallites, as indicated by the downward arrow. 

Our analyses compare collapse behavior within the frame- 
work of Fig. [T] We find that polymers and colloids behave 
similarly in many ways, and differences can be linked directly 
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Fig. 1 Schematic of potential energy versus temperature T for sys- 
tems with short-range attractive interactions cooled at varying quench 
rates T.— We compare the nonequilibrium and near-equilibrium 
crystallization dynamics of polymers and colloids by quenching 
systems at various rates and observing relaxation dynamics below 

^"cryst • 

to chain topology. Near-equilibrium crystallites possess sim- 
ilar static structure when T is scaled by the equilibrium crys- 
tallization temperature T cryst . Relaxation dynamics at fixed 
T/T cryst are qualitatively similar, but quantitatively quite dif- 
ferent for polymers and colloids. Using multiple measures of 
crystalline order and particle rearrangements, we quantify how 
restrictions on local motion imposed by the covalent backbone 
slow the approach to equilibrium by eliminating 'monomeric' 
relaxation mechanisms. Rearrangements in polymeric crystal- 
lites are required to be more cooperative, and hence are slower, 
than in their colloidal counterparts. We validate these results 
by examining whether these trends vary significantly with sys- 
tem size (i.e. chain length) TV. For N varying over a range 
typically studied in single chain polymer crystallization ex- 
periments, 12 increasing system size does not affect qualitative 
trends, but quantitatively strengthens them. 

The outline of the rest of the paper is as follows. In Section 
[21 we describe our model, simulation protocol, and metrics 
used to analyze our data. Section[3]presents results for thermal 
quenches over a wide range of \T\, nonequilibrium evolution 
of crystallite properties at fixed Tf < T cryst , and a detailed 
comparison of rearrangements in polymers and colloids. Fi- 
nally, in Section IU we summarize our findings and place them 
in context of other recent experimental and simulation studies. 

2 Methods 

Recent Monte Carlo simulations have employed advanced 
sampling techniques, such as topology-changing bridging 
moves, to investigate the equilibrium phase behavior 1,2 of 
single-chain polymers. In contrast, molecular dynamics (MD) 



simulations with physically realistic dynamics are better able 
to capture the complex, coordinated rearrangement events as- 
sociated with the glassy dynamics of crystallization. 

In our studies, both colloidal and polymeric systems con- 
sist of N identical spherical monomers that interact via the 
harmonic "sticky-sphere" potential shown in Fig.|2ja): 

k / r \ 2 ^ 

U havm (r) = I 6 2 ^ D ' T ° , (1) 

, r > r c 

where e is the intermonomer binding energy, D is the 
monomer diameter, and k = 1600e is the spring constant. The 
only difference between colloidal and polymeric interactions 
(inset to Fig.^a)) is that in polymeric systems, the monomers 
are linked into a linear chain connected by N — 1 permanent 
covalent bonds. 13 Different values for r c are used for cova- 
lently and noncovalently bonded monomers: r c c /D = oc and 
r™ c /D = 1 + ^/2e/fc, respectively. 

Newton's equations of motion are integrated using the 
velocity-Verlet method with a timestep dt = r/800, where 
m is the monomer mass and r = y / mD 2 /e. We deter- 
mined that this timestep was sufficiently small by examining 
^-dependence of the velocity autocorrelation function v ac (t) 
in simulations at high temperature; no statistically significant 
dependence was found for dt < r/600. Below, we express 
length scales, energies, times, rates and temperatures in units 
of D, e, r, t -1 , and e/fce, respectively. The temperature T 
is controlled via a Langevin thermostat with damping time 
10. A periodic cubic simulation cell with volume L 3 fixes the 
monomer density p. We present results for p = 0.01, which is 
in the dilute limit. 

Systems are initialized in random walk initial configura- 
tions (for polymers) and random nonoverlapping initial po- 
sitions (for colloidal systems). They are then equilibrated at 
high temperature, ksTi/e = 0.75, for times long compared 
to the time over which the self intermediate scattering func- 
tion decays to zero. At T = Ti, polymeric systems are in the 
good-solvent (self-avoiding coil) limit, and colloidal systems 
are in an ideal-gas-like state. Following equilibration, systems 
are thermally quenched at various rates T. In our ensemble- 
averaging approach, states from which thermal quenches are 
initiated are separated by times long compared to structural 
relaxation times, and thus are statistically independent. 

The quenches are either continued to T = or termi- 
nated at T = Tf < T cryst . In the latter case, we con- 
tinue the runs at fixed Tf. We choose Tf/T cryst = 7/8 
to suppress finite- N fluctuation effects associated with ther- 
mal broadening of the phase transition to crystallites, 15 i.e. 
Tf/T cryst < 1 — A^ -1 / 2 , yet allow sufficiently fast relaxation 
to be captured within the limits of available computational re- 
sources for the A^ considered here (40, 100, and 250). As we 
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will show below, this procedure yields particularly interest- 
ing results for nonequilibrium relaxation following quenches 
at moderate |T|. T//T cryst = 7/8 is also comparable to tem- 
peratures used in many experimental and simulation studies of 
crystallite nucleation and growth in supercooled colloidal and 
polymeric sy stems. -ISaii 

We will examine several order parameters to characterize 
nucleation, growth and rearrangements of crystallites as a 
function of temperature and time following thermal quenches 
to Tf/T cryst < 1. These order parameters are generated from 
the adjacency matrix A; Aij = Aji = 1 when particles i 
and j are in contact, i.e. when the position vectors satisfy 



<r r : 



, and otherwise. 



Our interaction potential (Eq. [TJ promotes contact- 
dominated crystallization dynamics. 18 For the large k/e and 
small r c employed here, the ground states for colloidal and 
polymeric systems 19 are simply the states that maximize the 
number of pair contacts N c = Ttj^Aij for a given iV.— ~— 
Further, our systems form crystallites possessing close-packed 
cores that increase in size as T decreases or equilibrium is ap- 
proached. 18 We therefore measure the number of close packed 
monomers N cp = H^S^jAij — 12) and degree distribu- 
tion P(r]), i.e. the probability for a particle to have r\ contacts. 
Note that Y%1 iP{i) = 2N C /N and P(12) = N cp /N. P{rj) 
contains more information than 7V C and N cp since it describes 
the high-77 cores of crystallites as well as their surfaces, where 
monomers naturally have lower 77 Q We will argue below that 
the combination of N c , N cp and P(rj) forms an effective set 
of "cry stal-agnostic" order parameters (in the spirit of Rein ten 
Wolde et.al.^). 

We will also present results for the adjacency matrix au- 
tocorrelation function Pamac(^ 5 t') as a function of t' after 
waiting various times t w following thermal quenches to Tf : 



^AMAC (tw,t') 



^ZiJ>i Aij(t w )Aij(t w + t ) 
^2i,j>i Aij(t w )Aij(t w ) 



(2) 



where the brackets indicate an ensemble average over inde- 
pendently prepared samples, and the total time elapsed after 
termination of the quench is t = t w + t' . For polymers, 
we exclude the contributions of covalent bonds to Pamac by 
summing over j > i + 1 rather than j > i in Eq. [2 Both 
^amac(^) an d the intermediate scattering function S(q,t') 
evaluated at qD ~ 2tt identify rearrangements of contacting 
neighbors, which control the slow relaxation processes in col- 
loidal systems with competing crystallization and glass tran- 
sitions.^^ To capture the glassy dynamics, we examine sys- 
tems with t w ranging over several orders of magnitude and 
t' > t w . 

*P(rj) is also closely related to the spectrum of eigenvalues of A 26 While 
these eigenspectra provide additional information on crystallite structure, we 
leave their examination for future studies of equilibrium systems. 



3 Results 

In this section, we compare the collapse and ordering dynam- 
ics of colloids and polymers using two protocols. Protocol 
(1) consists of decreasing the temperature from an initially 
high value Ti to zero using a wide range of thermal quench 
rates T. To analyze changes in structure with decreasing T, 
we measure the potential energy, N c , and N cp over the full- 
range of temperature for each quench rate. Protocol (2) con- 
sists of quenching systems from Ti to Tf = (7/8)T cryst us- 
ing a range of T and then monitoring structural order and 
rearrangement events within crystallites at T = Tf as sys- 
tems evolve toward equilibrium. We measure N c (t), N cp (t), 
P(rj] t) and PAMAc(t w ^ f ), to quantify evolution to more or- 
dered states, and also perform detailed studies of crystallite 
rearrangements in a "pre- terminal" relaxation regime where 
the systems slowly approach equilibrium. 

3.1 Protocol 1: Thermally Quench from High to Zero T 

Potential Energy: Figure [2] shows results for the scaled poten- 
tial energy U /Ne for colloidal and polymer systems quenched 
from T = Ti to zero. Panel (a) shows results at low \T\ for 
system sizes ranging over a factor of six in N, while panel (b) 
shows results for N = 100 over a range of \T\ spanning four 
orders of magnitude. All results are consistent with the gen- 
eral picture of Fig. [T] and illustrate both features common to 
polymers and colloids as well as differences arising from the 
presence of a covalent backbone. 

At the lowest \T\ considered (~ 10 -8 ), both colloids and 
polymers show sharp, first-order-like 1 transitions at corre- 
sponding T = T cryst . Because of the narrowness of the attrac- 
tive range of the potential well,^^ no intermediate liquid state 
(i.e., globules in the case of polymers) appears^^. In both 
cases, as in Fig.Q] the equilibrium transitions are from gas -like 
states to crystallites. No significant quench rate dependence is 
observable for T > T cryst , which indicates that all \T\ are 
sufficiently low to be near-equilibrium in this high tempera- 
ture regime. For \T\ ~ 10 -8 , polymers and colloids have the 
same energy at low T to within statistical noise, showing that 
this quench rate is slow enough to be in the near-equilibrium 
limit for polymers (i.e. \T*\ > 10 -8 for these systems.) The 
iV-dependence (panel (a)) shows only quantitative rather than 
qualitative differences. As N increases, values of U /Ne in the 
T — > limit decrease because larger crystallites possess more 
interior monomers. 

Compared to colloids, polymers have lower U /Ne for T > 
tcryst because of the permanent covalent bonds, which con- 
tribute ~ ksT/e — l. 30 They also have higher absolute values 
of T cryst (Table D. To zeroth order, T^jT* r f st ~ 3/2, a 
ratio which can be explained by a simple degree-of-freedom 
counting argument: while colloids have 3N — 6 nontrivial 
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Fig. 2 (a) System size dependence and (b) quench rate dependence of 
potential energy versus temperature from molecular dynamics simu- 
lations of polymers (solid curves) and colloids (dashed curves). Ver- 
tical dashed lines indicate T cryst . Panel (a) shows results for a single 
slow (|T| = 2.5 • 1(T 8 ) quench rate for N = 40 (red), N = 100 
(green), and N = 250 (blue). Panel (b) shows results for N = 100 
systems. Colors indicate \ f\ = 10~ 4 (red), 10~ 6 (green) and 10~ 8 
(blue). Results in (a) are averaged over 40 statistically independent 
samples, while results in (b) are averaged over 104 statistically inde- 
pendent samples for the higher two quench rates and 8 for the lowest. 
The inset to (a) illustrates the interaction potential (Eq. Q} employed 
in our simulations. 



degrees of freedom, the stiff covalent bonds in polymers act 
like holonomic constraints, reducing the effective dimension- 



ality of their phase space to 2 N — 5 so that 



pol 
cryst 



[(37V- 



6)/(27V - 5)]T c c r ° y ^ t ~ 3/2. 31 This observation helps motivate 
our (protocol 2) studies comparing relaxation of polymeric 
and colloidal crystallites at equal values of T/T cryst in terms 
of their different free energy landscapes. 

Table 1 Dependence of T crys t on N and topology. Our data are con- 
sistent with detailed analyses 8 predicting 0(N~ 1 ^ 3 ) finite-size cor- 



rections to T; 



pol 



cryst ' 



N 


rpcoll 

cryst 


rppol 

1 cryst 


rppol Irpcoll 

1 cryst / 1 cryst 


40 


0.237 


0.342 


1.44 


100 


0.240 


0.375 


1.56 


250 


0.239 


0.383 


1.60 



Below T cryst , dramatic quench rate dependence appears. 
Consistent with the general picture of Fig. [T] with increas- 
ing |T|, the T-dependent potential energy increases relative 
to the equilibrium-crystal limit. Figure |2b) shows results for 
N = 100. Similar trends are observed for the other N consid- 
ered, with the main difference being that |T*| increases and 
rate effects for \T\ > \T*\ become more dramatic with in- 
creasing TV. 

Contact number: We now discuss the temperature depen- 
dence of local measures of structure during quenches. Figure 
[3] shows N C (T) and N cp (T) plotted against T and T/T cryst 
for the same systems analyzed in Fig.[2tb). N c /N and N cp /N 
are particularly sensitive measures of crystallite equilibration. 
They increase monotonically with decreasing \T\ , but the vari- 
ations with rate are different for polymers and colloids. For 
fast quenches, polymers are more efficient crystal formers 
due to topological connectivity and the resultant cooperative 
dynamics. 7 An example of a cooperative polymer collapse 
mechanism not found in colloids is "chain pull-in"; a nucleus 
forms between chemically nearby monomers, and covalent 
bonds pull in segments that are chemically adjacent to the nu- 
cleus. 33 However, as exemplified by the \T\ = 10 -4 results, 
this mechanism produces only small ordered cores with highly 
disordered exteriors. Note that for the fastest quench rates, the 
collapsed structures possess N c /N > 3, which is the minimal 
number of contacts for mechanical stability 32 and corresponds 
to the isostaticity threshold shown in panel (a). 

In Fig. Eta), we show that |<97V C /<9T| at a notice- 
ably higher T for polymers than colloids, which indicates that 
large structural rearrangements cease, and polymers possess 
a higher T g i ass than colloids. 27 Polymer chain backbones fol- 
low tortuous paths through the interior of collapsed states^ 4 ^ 
and chain uncrossability suppresses large rearrangements , so a 
higher T g i ass is consistent with the expected slower relaxation 
dynamics for polymer crystallites. 
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Fig. 3 Measures of local structural order, (a) N c /N and (b) N cp /N, 
as a function of temperature during thermal quenches from T\ — 0.75 
to Tf = 0. (c) Number of close-packed spheres N cp /N plotted ver- 
sus temperature reduced by the crystallization temperature, T/T crys t. 
The line color, dashing scheme, and (N = 100) systems studied are 
the same as those in Fig. [2b). Vertical dashed lines indicate T cr yst 
for colloids and polymers, and the horizontal dashed line in panel (a) 
indicates the isostaticity 32 threshold (N c /N = 3). We define T g i aS s 
as the temperature at which rearrangements cease during a constant 
rate quench and the slope dN c /dT becomes close to zero. 



We find that measures of energy and contact number 
(U, N c , and N cp ) plotted against the scaled temperature 
T/T cryst for polymers and colloids all increasingly col- 
lapse as \T\ decreases. Figure [3c) shows N cp (T/T cryst ) 
for different quench rates. Over a broad range of T, the 
collapse of ensemble-averaged values of N C (T/T cryst ) and 
N C JT/T cryst ) for TV = 40 and N = 250 is of similar qual- 
ityQ This is a nontrivial result for the small crystallites consid- 
ered here, where ordering and surface-to- volume ratio changes 
with N and finite- size effects are in principle important. As 
shown in Table [2 terminal values of (N c /N) increase from 
~ 3.8 to ~ 4.7 as TV increases from 40 to 250 J3 and val- 
ues of (N cp /N) at T = increase even more (from ~ 0.15 
to ~ 0.35) as the volume-to-surface-area ratio of crystallites 
increases and they form larger close-packed cores. The col- 
lapse of U(T/T cryst ), 7V C (T/T cryst ), and N cp (T/T cryst ) de- 
spite these TV-dependent changes in order suggests that poly- 
mer and colloid crystallites, in equilibrium, occupy similar po- 
sitions on their respective energy landscapes at equal values of 
T/T cryst even though the absolute T are different. 

Table 2 Variation with N of (N c ) and (N cp ) during |T| = 2.5 • 10" 8 
quenches. Results are from the same systems depicted in Fig. 0a), 
and values of T crys t are given in Table Q] Results for (N c ) /N are 
consistent with leading order 0(N~ 1 ^ 3 ) corrections away from the 
N —> oo value (6). 



N 


T 


(N™ u /N) 


(NP° l /N) 


(N%"/N) 


(N?° L /N) 


40 


7 T 

g J- cryst 


3.73 


3.68 


0.143 


0.130 


40 





3.78 


3.76 


0.150 


0.146 


100 


7 T 

g J- cryst 


4.31 


4.24 


0.259 


0.239 


100 





4.33 


4.29 


0.265 


0.255 


250 


7 T 

g J- cryst 


4.67 


4.60 


0.362 


0.329 


250 





4.69 


4.65 


0.368 


0.343 



3.2 Protocol 2: Thermally quench from high T to Tf < 
T cryst and monitor relaxation at fixed Tf 

As shown in Table [2 polymeric and colloidal crystallites of 
the same N have nearly equal values of both (N c /N) and 
(N cp /N) at Tf = (7/8)T cryst . They also have similar struc- 
ture (i.e. P(rj)\ cf. Fig. [5]). This helps establish that compari- 
son of relaxation dynamics in polymeric and colloidal crystal- 
lites at fixed T = Tf is a reasonable approach to isolating the 



t Collapse worsens only slightly with increasing N, indicating that \T*\ in- 
creases slowly with N. Precise calculation of |T*(7V)| is outside the scope 
of this study, but more quantitative results for the TV-dependence of charac- 
teristic relaxation times is given in Section |3T2l 

X This variation is significant when one considers that the isostatic value 
(N c /N) = 3 can be attained by systems as small as TV = 16, 34 while 
the limiting value for infinite TV (corresponding to defect-free close-packed 
crystals) is {N c /N) = 6. 
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role played by topology in controlling the approach of crystal- 
lites to equilibrium. 

We now examine systems quenched at different rates and 
monitor their evolution as a function of time t following ter- 
mination of the quenches to Tf. Measures such as the evolu- 
tion of contact order and several measures of local and nonlo- 
cal rearrangements after different waiting times t w show that 
polymeric crystallites possess significantly slower long-time 
relaxation dynamics at equal values of T/T cryst despite their 
similar structure. Studies of system size dependence show that 
this "topological" slowdown is related to the more correlated 
rearrangements imposed by uncrossable covalent backbones, 
and strengthens with increasing N. 

Contact number: Fig. [4] shows the evolution of SN C , (N c ), 
and (N cp ) in polymeric and colloidal systems as a function 
of time t following thermal quenches to Tf at rates \T\ = 
10 -6 and 10 -7 . SN C is the root-mean- square fluctuation in 
the number of contacts averaged over an ensemble of collapse 
trajectories. The strong increase in (N c ) and (N cp ) and drop 
in 5N C after t ~ 10 5 for |T| = KT 6 in Fig. Ha) suggests 
that the crystal nucleation rate for colloids is r c « 10 -5 . In 
contrast, polymers do not show a rapid increase in the number 
of contacts (or concomitant decrease in the contact number 
fluctuations) at these thermal quench rates, showing that the 
crystal nucleation rate for (TV = 100) polymers is r p > 10 -5 . 

As shown in Fig. 0(b) and (c), at long times t > 10 6 both 
polymeric and colloidal systems show evidence of logarith- 
mic relaxation toward equilibrium. The logarithmic increase 
in (N c (t)) and (N cp (t)) is consistent with thermal activation 
over large energy barriers and transitions between metastable, 
globule-like states and near-equilibrium crystallites. 1,4 Relax- 
ation is also slowed by the increasing mechanical rigidity as- 
sociated with increasing 7V C .— The approach to equilibrium 
occurs through thermally-activated rearrangements of parti- 
cles associated with "soft modes" (cf. Figs. HH]), which are 
known to play a significant role in structural and stress relax- 
ation in supercooled liquids. 36 At larger t, the slopes (indi- 
cated by green solid lines) are clearly larger for colloids than 
for polymers; values of d (N c ) / d lnio t and d (N cp ) / d lnio t 
fit over the range 10 6,5 < t < 10 7 are given in Table [3] 

For t > 10 7 , values of d(N c (t)) /d\og 10 (t) and 
9{N cp (t)) / d\og 1Q (t) decrease for colloidal systems as they 
enter a terminal relaxation regime associated with ergodic ex- 
ploration of their full free energy landscape. The approach to 
the ergodic limit can be clearly seen in the vanishing of "his- 
tory" dependence for systems quenched at different rates, i.e. 
curves for different T overlap at large t. In contrast, for poly- 
mers, history dependence and faster logarithmic relaxation 
persist up to the maximum time t = 5 • 10 7 . This shows that 
polymers possess slower characteristic relaxation rates ,s s i ow 
at the same value of T/T cryst even though the crystallites are 
less mechanically rigid (since they have fewer contacts and are 
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Fig. 4 Nonequilibrium structural relaxation for polymeric (red) and 
colloidal (blue) systems after thermal quenches from T% > T crys t 



to T f 



(7/8)T c 



cryst 



at in 



10 



6 (heavier solid lines) and 10 7 



(lighter dashed lines). Results are averaged over 104 independent 
N = 100 samples and are plotted as a function of time t following 
the termination of the quenches: (a) Standard deviation in the number 
of contacts SN C = (jV c 2 ) - (N c ) 2 , (b) mean number of contacts 
(N c ), and (c) mean size of the close-packed cores (N cp ). In (b) and 
(c), the solid green lines show fits to logarithmic behavior at long 
times (cf. Table [5]); the lines are extended as a guide to the eye. 
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at higher absolute T). 

While observation of the crossover into this terminal re- 
laxation regime for polymers with N = 100 and TV = 250 
monomers is made unfeasible by the limitations of current 
computer power, in this paper we are primarily concerned with 
the nonequilibrium dynamics of crystallization and the loga- 
rithmic "pre- terminal" relaxation regime of crystallite growth. 
We now analyze the role of topology on relaxation dynamics 
within the pre- terminal regime by examining the evolution of 
more detailed measures of crystalline order. 




Fig. 5 Evolution of degree distributions P(rj). (a) P(rj) after 
\T\ — 10 -6 quenches to Tf. Curves show data averaged over times 
< t < 10 4 for colloids (dash-dotted) and polymers (dashed) and 
0.999 • 10 7 < t < 10 7 (colloids; light solid, polymers; heavy solid). 
Solid circles show results averaged over the 6 distinct N = 100 Bar- 
low packings. 37 Panel (b) shows AP(rj), obtained by subtracting the 
small-t data shown in panel (a) from large-t data for colloids (lighter 
blue line) and polymers (heavier red line). Inset: A common five-fold 
symmetric structure present on the surface of nonequilibrium crystal- 
lites. The green monomer has r] = 5. 

Degree distribution: Figure [5] shows results for the evolu- 
tion of P(rj) following \T\ — 10 -6 quenches. Polymeric and 
colloidal crystallites have similar degree distributions for in- 
termediate and high rj, indicating the crystallites' inner cores 
are similarly structured. However, covalent backbones pro- 
duce greater differences at crystallite surfaces. Polymer topol- 
ogy requires r\ > 2 for chemically interior monomers and 



r] > 1 for chain ends, while monomers in colloidal systems 
can have any degree of connectivity consistent with steric 
constraints (here, < r] < 12). Because of this difference 
in connectivity, polymer crystallites include a higher fraction 
of monomers with 2 < r] < 4 (panel (a)); this difference 
strengthens as t increases and systems approach equilibrium. 

Panel (b) shows AP(t?) = P(rj]t = 10 7 ) - P(rj]t = 0) 
to highlight the changes in P(r/; t) over an interval At = 10 7 
following the end of the quench to Tf . Notable features com- 
mon to colloidal and polymeric crystallites and depicted more 
clearly than in panel (a) include sharply negative AP(5) as- 
sociated with the annealing out of 5 -fold symmetric structures 
that tend to form on the surface of metastable crystallites (see 
inset), and sharply negative AP(ll) associated with anneal- 
ing out of defects within crystallite cores. The latter process is 
closely associated with the increase in (N cp (t)) shown in Fig. 
IU However, the key result shown here is that changes in P(rj) 
are uniformly smaller for polymers despite the higher absolute 
T. 

While the equilibrium P(rj) distribution is difficult to cal- 
culate for large N without resorting to advanced Monte Carlo 
techniques,^ one measure of evolving crystalline order in 
our systems can be obtained by comparing P(rj]t) to that 
of a known reference system. Barlow packings 38 are hard- 
sphere packings composed of layered hexagonal-close-packed 
planes; their three-dimensional order may be fee, hep, or 
mixed fec/hep, but they possess no defects (e.g. stack faults.) 
These are "reference" low energy structures for our model in 
the limit of steep hard-core repulsions and short-range attrac- 
tions. The solid circles in Fig. [3a) show P(rj) averaged over 
the six N = 100 Barlow packings. 37 

Our simulation data show that crystallites at the end of the 
pre-terminal relaxation regime possess Barlow-like order. In- 
terestingly, (N c ) is slightly higher (Fig. [4]) in the simulated 
systems than in the Barlow packings ((7V^ arlow ) = 421), 
while values of P(12) = (N cp /N) are similar. These effects 
are attributable to finite temperature, stiffness of core repul- 
sions and range of attractive interactions. Other differences 
between our systems and Barlow packings are attributable to 
small deviations from equilibrium, the fact that our crystallites 
may be stack-faulted, and (for polymers) entropic factors such 
as blocking. 21 That the above complexity can be understood 
by examining N c , N cp , and P(rj) shows that this combination 
of metrics constitutes an effective "crystal-agnostic" 16 mea- 
sure of orderJl 

In the above protocol 2 subsections, we have focused on 
results for TV = 100. System size effects are minor, e.g. 
slower logarithmic growth of crystalline order with increas- 
ing N and shift of P(rj) towards higher rj. Examining rear- 

§ Indeed, previous coarse-grained Monte Carlo studies of polymer crystalliza- 
tion^"— focusing on dense bulk systems have employed similar structural 
measures. 



This journal is ©The Royal Society of Chemistry [year] 



Journal Name, 20 1 0, [vol], 1 -\T3\ | 7 



rangements within crystallites provides additional insight into 
N- and topology-dependent effects on the glassy dynamics of 
crystallite formation and is discussed in the following subsec- 
tions. 

Adjacency matrix autocorrelation function: Next we exam- 
ine the decorrelation of contacts between neighboring parti- 
cles during the approach to equilibrium at Tf = (7/8)T cryst . 
We first examine effects of quench rate and waiting time on 
N = 100 systems, and then examine TV-dependence for evo- 
lution following slow quenches. The adjacency matrix auto- 
correlation function Pamac (t w > t') displays several important 
features (Fig. [6]): 

(i) : Polymeric rearrangement events are more frequent at 
low t' because of the higher absolute T. However, rearrange- 
ments are slower at large t' despite the higher absolute T. The 
slower decay arises from the covalent bonds in polymers that 
restrict the motion of monomer i to the plane tangent to the 
vector f^ + i — fi-i. Although contributions from permanent 
covalent bond contacts are excluded from the definition of 
^AMAcfc^')' m compact crystallites connectivity to chem- 
ically distant monomers produces long-range suppression of 
contact-breaking Q 

(ii) : Following the \T\ = 10 -6 thermal quenches, 
Pam ac (t w , t') for both polymers and colloids display strong 
t w -dependence as shown in Fig. Oa). We find an in- 
crease in the length and height of the low-t' "plateau" near 
Pamac (t w > t') ~ 0.95 with increasing t w , similar to the be- 
havior of the plateau in S(q, t w ,t') during the aging process 
in structural glasses. 27 Aging effects are stronger for colloids 
than for polymers because colloids are further from equilib- 
rium at the termination of the quenches (Fig.|4]a)). 

(iii) : For slower thermal quenches (|T| = 10 -7 ; Fig.[6](b)), 
similar but much weaker aging effects are present. Results 
for t w = and t w = 10 5,5 are indistinguishable to within 
statistical uncertainties for both polymers and colloids. Aging 
is delayed in part because the additional time to quench from 
Pcryst to Tf provided by the factor of 10 decrease in quench 
rate is larger than r" 1 ~ 10 5 and r" 1 > 10 5 (Fig.^a)), and 
in part because for t w > 10 6,5 systems have crossed into the 
preterminal (logarithmic) relaxation regime at t' =t — t w = 0. 
Both polymer and colloidal crystallites continue to slowly add 
contacts and close-packed monomers, and the ^-dependence 
should not vanish until equilibrium values of (N c ) and (N cp ) 
are reached. 

(iv) : At large t f , the adjacency matrix autocorrelation func- 
tion decays logarithmically. The crossover to logarithmic de- 
cay of Pamac (t w ■> t') corresponds to the pre-terminal regime 
of logarithmic growth of (N c ) and (N cp ) (Fig.^b) and (c)). 



This could be examined quantitatively by excluding successively more distant 
chemical neighbors (e.g. 2nd nearest, 3rd nearest), and by considering only 
chemically interior sections of polymers. A detailed analysis is left for future 
studies of near-equilibrium and equilibrium systems. 




10 2 10 3 10 4 10 5 10 6 10 7 



Fig. 6 The adjacency matrix autocorrelation function 
PAMAc(t w ,t') following (a) \f\ = 10" 6 and (b) 10" 7 ther- 
mal quenches. Protocols and systems for panels (a-b) are the same 
as in Fig. HI Results for polymers (colloids) are shown in red (blue). 
Pamac (t w ,t') is measured as a function of time t' after waiting 
times t w = 0, t w = 10 5 5 , and t w = 10 6 5 (solid, dashed, and 
dotted curves, respectively). The solid green lines in (b) show fits 
to logarithmic behavior at long times for t w = 10 6 5 and the black 
(+) sign denotes crossover into the terminal relaxation regime for 
colloids. Panel (c) shows how differences in Pamac (t w = 0;t ; ) 
between slowly quenched (|T| = 2.510 -8 ) polymeric and colloidal 
systems vary with N. Solid curves represent polymeric and dashed 
curves represent colloidal results, while colors indicate N = 40 
(red), N = 100 (green), and N = 250 (blue). Each curve in (c) 
represents an ensemble average over 40 statistically independent 
systems. 



For colloids and small t w , a decrease in the slope of this loga- 
rithmic decay corresponding to crossover into the terminal re- 
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laxation regime is indicated by the + symbol in Fig. Ob). No 
such change in slope occurs for polymers over the same range 
of t' . This is consistent with the idea that local relaxations in 
polymers are slower due to chain-connectivity constraints on 
rearrangements. 

(v): Figure Oc) illustrates the variation of Pamac(^w = 
0; t') with TV for slowly quenched systems. Characteristic con- 
tact decorrelation rates s cont decrease sharply with increasing 
TV; for example, the low-t' plateau lengthens with increasing 
TV, and the t' at which Pamac = 0.8 is 2-3 orders of mag- 
nitude greater for TV = 250 than for TV = 40. Since this 
TV-dependent decrease in s cont is similar for polymeric and 
colloidal crystallites (which, as we have shown, possess sim- 
ilar structure), we attribute it to the increasing contribution of 
crystallite cores where reneighboring dynamics are slow be- 
cause r\ is high and particles are more sterically constrained. 
Similarities and differences between polymeric and colloidal 
results are consistent with those expected from (i). For all TV, 
as in panels (a-b), polymers relax faster than colloids at low 
t' because Tf = (7/8)T cryst is higher, and slower at large t' 
because of topologically restricted rearrangement (cf. Figs. [7]- 
0) Both the "crossover" t' at which i^MAC = ^amac anc * 



the ratio Pamac I ^amac beyond this crossover time increase 
with increasing TV. 

In this paper we focus on the glassy dynamics of crystal- 
lite formation (where about half of the contacts existing at the 
termination of the quench have not been broken), not com- 
plete reorganization. Below, we show that there are significant 
differences between large-scale polymeric and colloidal rear- 
rangements in this regime. In the remainder of this section, 
we will focus on TV = 100 collapsed states generated using 
protocol 2 with thermal quench rate \T\ = 10 -7 . 

Statistical comparison of rearrangements in polymeric and 
colloidal crystallites: We describe rearrangement events us- 
ing the two-dimensional parameter space (TV c ,TV cp ), where 
tangent- sticky- sphere polymers and colloids have the same in- 
herent structures.— In Fig. [7] we show the probability distri- 
bution P(ATV C , AN cp ) for crystallite rearrangements to add 
ATV C contacts and ATV cp close-packed particles in crystallites 
over time intervals At = 10 3 Ui| Results are calculated for the 
range 10 6 ' 5 < t < 10 T where both colloids and polymers are 
in the preterminal relaxation regime. P(ATV, 
portional to the integrated rate 



c ATV cp ) is pro- 



i?(ATV c ,ATV cp ) 



for all transitions that add ATV r contacts and ATV 



cp 



cp 

(3) 
close- 



packed particles in crystallites originally posessing N® con- 

| Note that for the physically reasonable values m = 10 _24 kg, D = lnm, 
and e ~ 10k bT at room temperature, r ~ 5ps, and timescales ~ At can be 
probed by neutron spin echo experiments, e.g. for the purpose of characteriz- 
ing protein dynamics.— 



tacts and N® close-packed particles, i.e. P(ATV C , AN cp ) ~ 
P(ATV C , AN cp ) / At. This becomes exact in the limit 
R(AN c ,AN cp )At 0. Thus, the data in Fig. [7] provides 
a basis for comparing free energy barriers and transition rates 
in these systems. 

Figure [7] and Table [3j which characterize the shape of 
the distribution P(ATV C , ATV cp ), illustrate the dramatic differ- 
ences between polymeric and colloidal rearrangements during 
logarithmic relaxation. Polymeric rearrangements are signif- 
icantly more correlated than those for colloids. For example, 
the correlation coefficient 



( (ATV C - (ATV C )) (ATV cp - (ATV cp )) ) 
( (ATV C - (ATV C )) 2 (ATV cp - (ATV cp )) 2 ) 



(4) 



where the averages are taken over all rearrangements, is larger 
by a factor of 2.5 and the cross correlation (ATV c ATV cp ) 
is larger by a factor of 7.4 for polymers compared to col- 
loids. The second moments of P(ATV C , ATV cp ), ((ATV C ) 2 ) 
and ((ATV cp ) 2 ), are also larger for polymers than colloids. 

Table 3 Statistical analysis of the data presented in Fig. [7] The top 
two rows are calculated by fitting to data in Fig. |4]in the pretermi- 
nal regime (10 6 5 < t < 10 7 ) The middle column shows data for 
colloidal rearrangements excluding "floaters". 



Quantity 



Colloids Colloids (NF) Polymers 



d(N c ) /d\og 10 t 

d{N cp )/d\og 10 t 

<ATV C 2 ) 

(ATV c 2 p > 

(ATV c ATV cp ) 

c 



6.3 


6.3 


1.4 


2.1 


2.1 


0.7 


9.1 


2.2 


23 


0.41 


0.32 


3.2 


0.66 


0.39 


4.9 


0.23 


0.34 


0.58 



An interesting feature of P(ATV C , AN cp ) is that it is 
nonzero in quadrants II and IV (i. e. negative AN c AN cp ). 
Rearrangements occur in which the number of close-packed 
monomers decreases but the number of contacts increases 
overall, and vice versa. These are less likely for polymers than 
for colloids since covalent bonds impart greater cooperativity 
to rearrangements in polymers. P(ATV C , AN cp ) is reasonably 
well fit by a 2D Gaussian functional form 



P*(ATV C , ATV cp ) OC e -E{AN c )^FAN c AN cp -G(AN c: 



(5) 



where E is similar for polymers and colloids, F is larger for 
polymers, G is larger for colloids, and —Fj (2EG) is larger 
for polymers. These results are consistent with our finding 
that c is larger for polymers than for colloids. 

The fit P*(ATV C , ATV cp ) does, however, fail to cap- 
ture several key features of rearrangements. In particular, 
P(ATV C , AN cp ) has a sharper peak at the origin and "fat" 
(wider than Gaussian) tails. The fat tails of P(ATV C , AN cp ) 
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Fig. 7 Probability distributions P(AN C , AN cp ) for changes in the number of contacts and number of close-packed monomers over time 
intervals of At = 10 3 following \T\ — 10 -7 thermal quenches from > T cryst to T//T cryst = 7/8 for (a) colloids, (b) colloids with 
"floaters" excluded, and (c) polymers. The different colored regions indicate bins in probability that differ by factors of 10. Results are for the 
same systems analyzed in Figs. |4] and Ob). To capture the "pre -terminal" logarithmic relaxation regime, only results for 10 6 5 < t < 10 7 are 
presented. 



are dominated by events where a monomer escapes from or 
rejoins the crystallite (for colloids), or a chain end becomes 
"floppy" (for polymers). These events tend to be associated 
with rearrangements with large AN C and small AN cp for col- 
loids, and large A7V C and AN cp for polymers (the difference 
being that covalent bonds impart a large degree of coupling 
between the exterior and close-packed interior of crystallites). 
While rare, these events may dominate large rearrangements 
in which many changes of contacts occur. 

To quantify this effect, we analyze colloidal rearrangements 
with "floaters" excluded, i.e. where all particles possess at 
least one interparticle contact both before and after the rear- 
rangement. Since monomers in polymers always have at least 
one contact, excluding floaters provides a more direct method 
to quantify the role of topology in controlling the rearrange- 
ments occurring during slow crystallite growth. Figure 0b) 
and the "Colloids (NF)" data in Table [3] present this analysis. 
The "no-floater" subset of colloidal rearrangements possesses 
significantly smaller fluctuations (in terms of (AAT?) and 
(AiV c 2 p )) and is more correlated (in terms of (AN c AN cp )). 
However, no-floater colloidal rearrangements remain signifi- 
cantly less correlated than those for polymers. Thus we claim 
that the distributions P(A7V C , AN cp ) differ dramatically for 
polymers and colloids (i.e. are more correlated for polymers) 
precisely because their free energy landscapes are different, 
which in turn is a direct consequence of the presence or lack 
of covalent bonds. 

The large width of the distribution P(A7V C , AN cp ) implies 
that the dynamics of crystallite growth in the pre-terminal re- 



laxation regime are highly heterogeneous, especially for poly- 
mers. One source of dynamical heterogeneity, as suggested 
by Fig. [71 is surface effects. Rearrangements at the surfaces of 
crystallites (where monomers possess low rj) are qualitatively 
different than those occurring within the close-packed cores. 
To illustrate the role of surface effects and covalent bonds, we 
now compare typical rearrangement events characteristic of 
those occurring in (N = 100) polymeric and colloidal crystal- 
lites. 

Visualization of typical rearrangements in polymeric and 
colloidal crystallites: Figure [8] shows a colloidal rearrange- 
ment with AN C = 8, AN cp = 2 ((a) and (b)) and a polymeric 
rearrangement with AN C = 7, AN cp = 4 ((c) and (d)). These 
are roughly equiprobable, with log 10 P(AN C , AN cp ) « 
-3.5, and the ratio AN cp /AN c w c. (See Table EH) Be- 
fore the colloidal rearrangement, the red monomer has only 
one contact; during the rearrangement, it settles into a groove 
on the surface of the crystallite and adds three new contacts. 
At the same time, the right hand side of the crystallite under- 
goes a large rearrangement, which improves its stacking order, 
and five other contacts and two other close-packed atoms are 
added elsewhere. It is likely that the initial "looseness" of the 
red monomer destabilizes the crystallite and gives rise to the 
soft mode that causes this large rearrangement. 

For polymers, the picture is qualitatively similar, but is fur- 
ther complicated by the fact that the interior and exterior of 
crystallites are topologically connected. Large rearrangements 
are often initiated when a chain end is at the surface and rel- 
atively loose. Note that in Fig. [8] (c) and (d) the color varies 
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AN C = 8, AN C p = 2 



"Loose" particle 
with 1 contact 
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adds 3 contacts 
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Polymers 

AN C = 7, AN C p = 4 



Free volume pocket 
Before 




Pocket collapses, creates 
new close-packed atoms 



After 




Fig. 8 Visualizations of typical moderately- sized rearrangement 
events for polymeric and colloidal crystallites. Panels (a) and (b) 
show the colloidal crystallite before and after a colloidal rearrange- 
ment event with AiV c = 8 and AN cp = 2 and (c) and (d) show 
the polymeric crystallite before and after a polymeric rearrangement 
event with AN C — 7 and AN cp = 4. The orange color shading in 
(c) and line in (d) show where new close-packed particles are added. 



from red to green to blue as the monomer index increases from 
1 to 100. Before the rearrangement, the chain end including 
monomer 100 is relatively loose (with only 3 contacts), and a 
pocket exists in the crystallite with above average free volume. 
During the rearrangement a segment of the polymer including 
this end executes a "flip" that collapses the pocket. The chain 
end monomer adds three contacts, and four other contacts and 
close-packed atoms are added elsewhere. 

Covalent bonds suppress large rearrangements less when 
chain ends exist on the exterior of crystallites. In the rear- 
rangement event depicted in Fig. [8] (c) and (d), the path of 



the covalent backbone through the crystallite is not particu- 
larly tortuous — it proceeds in a relatively orderly fashion from 
upper-right to lower-left. Entropic factors such as "blocking" 
suppress the probability for chain ends to exist in the inte- 
r ^ or ^39 othe^ise large rearrangements would be even fur- 
ther suppressed. 

In summary, despite qualitative similarities, polymer topol- 
ogy produces the large quantitative differences in slow crys- 
tallite growth and rearrangements illustrated above in Figs. 
HH] and Table [3] Even larger and rarer polymeric rearrange- 
ments than those depicted in Fig. [5] involve cooperative re- 
arrangements of sub-chains that do not include chain ends. 
Such large-scale rearrangements occur within the interior of 
the crystallite; their initiation requires a large (negative) AN C 
and AN cp , and hence they possess large activation energies. 
These are the slowest relaxation mechanisms, and it is likely 
that they control the approach to the ergodic limit in polymer 
crystallites. 

4 Discussion and Conclusions 

In this article, we compared the crystallization dynamics of 
single-chain polymers and colloids. The use of model sys- 
tems with hard-core-like repulsive and short-range attractive 
interparticle potentials yielded contact-dominated crystalliza- 
tion and allowed us to characterize crystalline order via mea- 
sures such as the number of contacts N c , the number of 
close-packed particles N cp , and the contact degree distribution 
P(rj). Our use of a model in which covalent and noncovalent 
bonds have the same equilibrium bond length yielded the same 
low energy structures for polymeric and colloidal systems, al- 
lowing us to isolate the role of chain topology on the dynamics 
of crystallite formation and growth. 

Particular attention was paid to the effect of thermal quench 
rate. Slow thermal quench rates yield first-order like transi- 
tions to crystallites at T = T cryst . The ratio of T cryst for poly- 
mers to that for colloids can be obtained roughly by counting 
degrees of freedom. Comparison of polymeric and colloidal 
crystals at equal values of T/T cryst showed they possess simi- 
lar structure (i.e. N c , N cp , and P(rj)), and thus occupy similar 
positions on their respective free energy landscapes, despite 
significantly different absolute T. Higher quench rates yielded 
rate-dependent effects and glassy relaxation from partially dis- 
ordered to more ordered configurations. While the marked 
slowdown in dynamics at T = T cryst and consequent rate- 
dependent glassy behavior for crystal-forming systems pos- 
sessing phase diagrams like Fig. |2a) is understood in terms 
of a crossover to potential energy landscape dominated dy- 
namics 44 with decreasing temperature, the role of covalent 
backbones (and consequently, different energy landcapes) in 
producing the strongly differing nonequilibrium responses for 
polymers and colloids reported in this paper has not been pre- 
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viously isolated. We showed that although polymer crystal- 
lites nucleate faster because of the cooperative dynamics im- 
parted by their covalent backbones, chain connectivity slows 
their relaxation towards maximally ordered structures. 

Crystallites can rearrange in many different ways (e.g. with 
different changes AN C and AN cp in the number of contacts 
and close-packed monomers, respectively). By measuring the 
transition probability P(A7V C , AN cp ) in the regime where the 
degree of crystalline order exhibits slow logarithmic growth, 
we characterized how the rare collective rearrangement events 
which control the slow approach of crystallites to equilibrium 
are affected by polymer topology. Significant differences be- 
tween P(AN CJ AN cp ) in polymeric and colloidal crystallite 
formation are attributable to the increased cooperativity of re- 
arrangements required by the covalent backbone. 

Strong finite size effects have been observed in both equilib- 
rium and nonequilibrium polymer-collapse studies. 217 45 Here 
we examined system size effects using measures of order such 
as N c and N cp that vary rapidly with the number of particles N 
due to concomitant variation in the ratio of crystallite surface 
area to volume. This variation did not change any of the qual- 
itative features reported above, and quantitative differences 
were as expected; dynamical slowdown of relaxation associ- 
ated with restricted motion imposed by the covalent backbone 
strengthens with increasing N. 

It is well known that bond-angle interactions play a signifi- 
cant role in controlling crystallization of most synthetic poly- 
mers. While this study considered flexible chains, it serves as 
a basis for future studies of more realistic models by elucidat- 
ing the role polymer topology plays in controlling the glassy 
dynamics of crystallization. Our results may also be directly 
applicable to understanding the collapse behavior of flexible 
"colloidal polymers", which have recently been shown to self- 
assemble into tunable, compact nanostructures,— as well as 
very flexible natural polymers such as single stranded DNA. 46 
Future work will examine how crystallization and packing are 
affected by semiflexibility, as well as effects of topology on 
the dynamics of equilibrium crystallites. 
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